Regular collision of dilatonic inflating branes 
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We demonstrate that a two brane system with a bulk scalar field driving power-law inflation on 
the branes has an instability in the radion. We solve for the resulting trajectory of the brane, and 
find that the instabihty can lead to collision. Brane quantities such as the scale factor are shown 
to be regular at this collision. In addition we describe the system using a low energy expansion. 
The low energy expansion accurately reproduces the known exact solution, but also identifies an 
alternative solution for the bulk metric and brane trajectory. 



I. INTRODUCTION 

In the search for a theory that can describe 
both quantum mechanics and gravity, string theory 
(which describes matter as vibrating strings moving 
in 11 dimensions) seems a strong candidate. It has 
been the subject of much research, and many differ- 
ent realisations of string theory have been found to 
be possible. This being so, experimental constraints 
on the theory are vital. It is hoped that the next 
generation of particle accelerators now being built 
will reach energies where some sign of extra dimen- 
sions may be seen, while laboratory gravity exper- 
iments at low energies can put constraints on the 
length scale of any extra dimensions. In cosmology, 
we can probe both the highest energies, through un- 
derstanding of the early universe, and the largest 
scales available to us; thus an understanding of how 
string theory and extra dimensions affects cosmology 
is very valuable. 

The fact that we observe only three spatial di- 
mensions has to be explained away in string the- 
ory - this can be done by compactifying the ex- 
tra dimensions to very small scales, or by attach- 
ing all the strings corresponding to normal matter 
to a four-dimensional surface called a brane. Five- 
dimensional models based on this latter method of 
dimensional reduction were presented in two papers 
by Randall and Sundrum in 1999 0. Much work 
has been done on cosmological generalisations of 
these models since then. Early universe inflation 
in braneworlds has been extensively studied, mostly 
driven by some kind of inflaton field, either localised 
on the brane with the other matter fields I3| , or liv- 
ing in the whole bulk 3, 4] . In this paper we will be 
considering a model of bulk-driven inflation. 

In Isl] , Koyama and Takahashi found an exact so- 
lution for a five-dimensional bulk metric where a 
bulk scalar field with non-BPS exponential potential 
drives power-law inflation on a single brane. This 
scenario is particularly attractive as the cosmologi- 
cal perturbations, normally an intractable problem 
in braneworld scenarios, can be solved exactly (Qj. 

In 0, Mukohyama and Coley developed the sce- 



nario to include two branes. They found that for any 
value of the brane tension on the second brane there 
is one location where the second brane will remain at 
a constant bulk coordinate away from the first brane. 
This is somewhat analogous to the scenario with two 
de Sitter (dS) branes and an empty Anti-de Sitter 
(AdS) bulk studied by Gen and Sasaki H. However 
unlike the static de Sitter model, the presence of 
the scalar field means that the constant coordinate 
separation does not lead to a static equilibrium as 
the proper distance between the two branes is time 
dependent. 



The stability of this two-brane model against a 
wide class of metric perturbations was shown in 0. 
Although it was shown that there were no unstable 
bulk metric perturbations, including perturbations 
of the bulk metric introduced by radion fluctuations, 
the radion itself was not calculated. Thus a scenario 
with a brane moving in the fixed bulk space-time 
was not excluded by the analysis of Such a sce- 
nario would be the natural analogue of the radion 
instability found in the two-dS brane model in j^. 
In Section FlI Bl we will examine the case where the 
second brane is allowed to move in the fixed bulk 
metric, and solve the junction conditions giving its 
trajectory. We find that the second brane position 
is indeed unstable, with any displacement leading to 
the second brane going to infinity or colliding with 
the reference brane. However, it is shown that all 
brane quantities are completely regular at this col- 
lision. This is in contrast with the results of Web- 
ster and Davis |^, who examined similar non-BPS 
braneworld models with more general scalar field po- 
tentials, and found that the scalar field typically be- 
comes divergent at a collision. 



The stability of this scenario to homogenous met- 
ric perturbations is examined in a low energy expan- 
sion in Section These results reinforce the sta- 
bility found in 0. We also demonstrate the ability 
of the low energy expansion to reproduce a known 
exact solution in the correct limit. 
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II. EXACT SOLUTION 
A. Previous Results 

The scalar field $ in the bulk has potential j^: 
A($) =n-^(-+s) a2e-2^'"'*(--*) , (1) 



and induced potentials on the two branes: 

Aoa($) = V2K-Vo,ie-^^'^*(^'*), (2) 

where $ is taken to have a separable form: 

$(z,i) = (/.(t) + S(z) . (3) 

The parameter A is defined in terms of the coupling 
constant b as 



46^ 



(4) 



Then the part of the bulk potential proportional to 
A is the BPS potential, and tuned to the brane po- 
tentials such that the metric induced on the branes 
will be static (this is equivalent to the Randall- 
Sundrum tuning in the model with AdS bulk). The 
parameter 5 is a free parameter giving the devia- 
tion from this BPS tuning, and is what allows the 
time-dependent evolution on the brane. 
The action for this model is given by 



We assume that —8/3 < A < —2 so that all expo- 
nents are finite, and that A/8 + S < Oto ensure A < 
0. Then we note that for (5 7^ a brane at position 
z = zo, determined by 7i(zo) = 1, undergoes power 
law inflation with scale factor a{t) = {Hot)'^/''^'^+^^; 
and for a brane at this location bulk time t coincides 
with brane proper time. In the limit S —^ we have 
that H = Ho = 0, giving a scale factor of 0, demon- 
strating that the power law inflation is driven by the 
scalar field and by the 5 part of the scalar field in 
particular. 

A second brane may be introduced at arbitrary 
constant 2; = zi as in Q. Proper distance between 
the branes is then given by 



H, 



ot I 

J Zo 



(13) 



In the case of two dS branes the equilibrium position 
has constant proper distance; the time-dependence 
of the power-law inflation on the brane causes a 
simple time-dependence for the equilibrium proper 
distance. As Hq is positive, once the branes have 
emerged from the curvature singularity at t = they 
are always moving apart. 

The junction conditions at the first brane are au- 
tomatically satisfied; at the second brane they give 
a single equation relating zi and the two brane ten- 
sions, (To and (Ti: 



S = 



(5) 



Taking the bulk metric ansatz: 

the Einstein equations are solved by 0: 



(6) 





(7) 




(8) 




(9) 




(10) 



with 
Ho 



3A 



3(A + 2) 
and for A < 0, H(z) is given by: 



H, H = -{A+2)\l~-ao (11) 



niz) = \/-l - ^sinhiJz 



(12) 



In proper time on the second brane the scale factor 
gives power law inflation with the same exponent 
but a different constant of proportionality. 



B. Perturbing the equilibrium position 

We now relax the assumption that the second 
brane has to be at constant z and perturb the brane 
from its equilibrium position found above. The gov- 
erning equations are derived from the junction con- 
ditions and prove not to be analytically soluble, 
though a numeric solution can be obtained. Nu- 
meric solutions shown here have initial conditions 
close to, but slightly perturbed from the constant z 
brane position. 



1. Junction Conditions 

The power-law inflation demonstrated in Section 
III Al ls induced by the action of the bulk scalar field 
on the brane, with the brane being empty apart from 
the scalar field potential Q . We now allow for mat- 
ter with energy density p and pressure p on the sec- 
ond brane. Then we find the junction conditions and 
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the scalar field matching conditions (given generally 
in 0) at the second brane become: 



1, \21-3/2 ^^w 



W 



-W'Hotzl + 2iJoii + — (15) 



Hut 



6 



■{Ai($)-2p-3p} 



= -y{Ai($)+p}, 

[l-Wz,f]-"\-'^ 

, 2HqZi 



2HoZi 



W 



Hot (3A 



(16) 



(3A + 8) 
W ■ 



(17) 



We can see that (dSl) and ^ are inconsistent 
unless p = 0. The bulk metric ansatz and the form 
of the scalar field matching condition we have cho- 
sen do not allow for matter on the brane. It would 
be possible to add matter by taking a more general 
metric ansatz, or alternatively allowing a coupling 
of the scalar field to matter to change the form of 
equation El (a-s discussed in jl^). Here, however, 
we will investigate the case with no matter on the 
branes. 

When p is taken to be zero, equations (|16|l and 
(|17|l are equivalent, and equation H15|l is equivalent 
to the derivative of these two. 

To evaluate the effect of any deviation from the 
zi = const scenario previously considered we make 
the perturbation zi 'zi-\-5z, where zi is that value 
of zi that satisfies equation (|14|l and 5z is small. 
Then the linearised expansion of equation (|16|l gives: 



—5{e^W') 

Hot 



3A + 8 ^ ~ 6^^^ H^t 



-6z. 
^ (18) 

when we substitute the solution for W and ^ given 
(in the case A < 0) by equations - l(T^. this 
simplifies to 



3(A + 2) Sz 
'(3A + 8)T 



= Sz. 



(19) 



This has solution 5z cx t-3(A+2)/(3A+8)^ p^j. ^-^e 
range of A we are considering, the exponent of t 
is positive, and hence the radion is unstable. We 
note that this solution agrees with the inverse fourier 
transform of equation (49) of reference in the case 
where all bulk metric perturbations are zero. 



2. A numeric solution for the brane trajectory 

To go beyond a linear analysis we solve equation 
(O for ii: 



zi = 



72(3A + 1 



n 



72H^n^ + afiSA + t 
V2 



3(A + 2) 



(20) 




(7f(3A 



72H^W 



(3A + 8)2 

9(A + 2)2 H^H^ 



We introduce coordinates normalised with respect 
to reference brane tension ao- 



H 



H 



11 



t aot , z 



aoz. 



In order to obtain a numeric solution of the differ- 
ential equation H2U|) we choose values for the param- 
eters as: 



A 



4 ' 4 



H = 



1 

12 



Ho^—,ai^-2. 
3d 



(21) 



These have been chosen to be consistent with the 
definitions and inequalities presented in Section fll Al 
Then the position of the reference brane is given by 



zo = 12sinh-^ (2V2) ~ 21.2, 



(22) 



and the tuned second brane position giving constant 
zi is (from equation H14|l ) 



zi = 12cosh"^ (6) ~ 29.7. 



(23) 



From equations H12|l and H2U|) we get the equation 
for zi{t): 

_ 36 [sinh(zi/12) cosh (zi/12) - 6^35] 
^'~T [{cosh(zi/12)}2 + 35]i ■ 

Figure ^ shows numeric solutions of this equa- 
tion with initial conditions chosen so that the three 
curves show the trajectory of the second brane when 
it is initially at or slightly to either side of, the tuned 
position calculated above in equation (I23() . The ref- 
erence brane position is also shown, and one can see 
that when the second brane is perturbed towards the 
reference brane, they will eventually collide. 



3. Scale factor on the moving brane 

As shown in [Tlj | the five-dimensional Einstein 
equations can be projected onto a brane to obtain 
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four-dimensional effective equations. In the case of 
a bulk scalar field with no matter on the branes, 
the effective four-dimensional equations on the brane 
are [H: 

^'^G,, = J-ale-^^^'-^ + \>^'T,^m - E,. , (25) 

where T^^{<^) = D ^'^D^'^ - \q^,^(D'^f , the covari- 
ant derivative on the brane is and the induced 
metric on the brane is q^^. E^^, is the projection of 
the bulk Weyl tensor onto the brane. 

The Weyl part of these effective Einstein equa- 
tions can be broken up into a piece dependent on $ 
and a piece independent of it as follows [T^ : 



E„ 



1 



(26) 



E^ii, is necessarily traceless, giving 



where Sv^&'t-Fi is the trace of the $ part of -E^jy, 
given by 



(28) 



35 



20 



15DD 
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FIG. 1: Second brane position zxify plotted against bulk 
coordinate time t. Initial conditions are 2i(0) = z\ 
(dashed line) and 2:1(0) = 21 ±0.1 (dotted lines). The ref- 
erence brane position is also shown (solid line). The so- 
lution for a brane perturbed towards the reference brane 
leads to a collision at time tc — 1920, or brane proper 
time Tc ~ 919. 



The bulk metric considered here is further con- 
strained to have F$ = = F^y, where = gives 
the wave equation for (/). 

The symmetries of the bulk metric combined with 
the wave equation for the scalar field determine the 
induced metric and scalar field on any brane in the 
bulk metric. There are insufficient degrees of free- 
dom to allow for several independent solutions, with 
the result that these quantities will have the same 
r-dependence on a moving brane as on a fixed one. 
This result is demonstrated below using our numeric 
solution for the brane trajectory. 

From the differential equation (|24() and the form of 
the bulk metric we obtain the induced metric along 
the brane trajectory as: 



= -dr^ + 256[sinh(zi/12)]-i6/3 



16/5 



dx , 
(29) 



where proper time r along the brane trajectory is 
related to bulk coordinate time t by 



dT 



16 



6 cosh (01/12) + V35sinh (zi/12) 



(27) dt [sinh(zi/12)]8/3 



[{cosh(zi/12)}2 + 35] 
(30) 

Choosing the numeric solution for zi with the second 
brane moving towards the reference brane, equation 
(|30() can be solved numerically and inverted to get 
t{r). 

From the induced metric (|29|l we have that the 
scale factor ai on the second brane is 



ai 



16 



5t 



[sinh(zi/12)]8/3 V36 
and, using (|30|) . we have 



8/5 



(31) 



1 dai 

Hi = — = 

ai dT 



lot 



sinh(zi/12)]'5/3 (32) 



A plot of ai (r) is shown in Figure|21 A calculation 
of the quantity tHi (r) in the numeric solution shows 
that it is constant at a value of 2/(3A -|- 8) = 8/5, 
demonstrating that the scale-factor is power-law 
with the same exponent as the non-moving brane. 



C. Constraints on the bulk metric 

In the above we have used a particular constrained 
form of the bulk metric whose symmetry properties 
discounting any possibility that the brane cosmol- 
ogy or the movement of the brane through the bulk 
affect the form of the bulk metric. In particular, the 
bulk scalar field and the bulk Weyl tensor are con- 
strained to symmetric forms. This symmetry may 
be described in terms of the two quantities F^ and 
defined above. 
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distance being given by the formula 

L = e^'"^d(i) . (33) 

In this section we wih use the variables and pa- 
rameters of JJj] , but to ensure we are describing the 
physical system set up previously, the parameters 
must be given by a = ao, = -25a'^e-'^'^^'t' , 

U{4>) = and U{4>) = d'y^/K^e-^*"^ where a' = 
a^) — (Ji\ then the action given in [T^ reduces to that 
of equation (O . 

The low energy expansion given in 17] relies on 
the assumption 

a > K^p, K^p, a~^V, K^U, K^ij , (34) 

which in our case implies 



FIG. 2: Scale factor ai on the second brane plotted 
against brane proper time r for the trajectory perturbed 
towards the reference brane. We can see that there is 
no singularity or other feature around the collision point 
Tc ~ 919. 



(5 < 1, ct' < cr. 



(35) 



In the d = dif — constant, zi — constant case we 
can derive the coordinate transformation between 
the exact metric given by equations (|7|)- ((Tn)l and 
the approximate one of |17| as: a; — > x , t —> t and 



Specifying the bulk metric ensures that this sym- 
metry will not be disturbed if the brane is perturbed 
and starts to move. This results in a first order evo- 
lution equation for the radion, where other works 
have found second order equations 14]. Given an 
initial value for the radion, the brane velocity is fixed 
by the evolution equation. To set arbitrary initial 
position and velocity, we must allow for some per- 
turbation of the bulk metric parameters. 

Solving the full system of Einstein equations for 
a moving brane would be far too complicated to be 
practical, but if we use a low energy expansion we 
can examine the possibility for a more general bulk 
in the low energy regime. 



12 



V2{3A + 8)d 



(36) 



A 



(3A + 8) 
3(A + 2) 



sinh (Hz) 



(valid at low energies and small S only). From this 
transformation and the low energy conditions H35|l . 
we can derive the condition on the exact solution 
that Hz < Hzi ^ 1. We suppose that this condi- 
tion would also hold when d is not constant. 



A. Linear approximation 



III. LOW ENERGY EXPANSION 

The 4-dimensional effective theory valid at low en- 
ergies for a system with 2 branes is derived in Ts','!^ 
and applied to a bulk scalar field in |0| . We will sub- 
stitute the potentials from equations ^ and into 
this low energy formalism to determine the validity 
of allowing free brane motion while constraining the 
bulk to the static metric given in 

In this low energy theory, there are essentially 
three moduli fields d, 4> and d. a and 4> correspond 
to the scale-factor and the scalar field on the refer- 
ence brane, and d encodes an arbitrary variation of 
the proper distance between the branes; the proper 



The effective equations on the brane are given by 
equations (58)-(60) in . To evaluate the effect 
of some small deviation from the zi = constant, 
d — constant scenario we make a linear perturba- 
tion d = d., -\- 5d where 5d is small compared to d, . 
We assume that a' is the tuned brane tension to keep 
d at constant d*: 



4(5cr 



(1 



\(6A+12)/(3A+8) 



(37) 



and that any deviation of (p and a from their values 
in the exact solution is also of order 5d or smaller. 
Under these assumptions, the evolution equation for 
d decouples from the governing equations for (p and 
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d, and gives a second order ordinary differential 
equation, 



3A + 14\ 5d 



3A 



t 



9A(A + 2) 



(3A + i 



Sri 

^ = 0. (38) 



Substituting — into the evolution equation for (j) 
obtained in the low energy expansion and combining 
with the radion evolution equation H38|l gives a first 
order equation 



Sri. 



3(A + 2) 5d 
"(3A + 8) t 



(39) 



The second order evolution equation H38() has two 
independent solutions: 



Sri OC t-3(A+2)/(3A+8) ^ 
Sri OC t3A/(3A+8)^ 



(40) 
(41) 



The growing mode H4U|I satisfies the F,j, — condi- 
tion H39() . while the decaying mode (|41|l does not. If 
one chooses initial conditions for the brane motion 
satisfying = 0, they specify solely the growing 
mode, so a solution that starts with F,p = will re- 
main that way for all time. The one-parameter fam- 
ily of solutions Section lll BI must lie along this grow- 
ing mode trajectory. Any solution that also contains 
some decaying mode solution will not have F^ = 0, 
but it will tend towards a trajectory that does, im- 
plying that Frf, = is in some sense an attractor. 
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FIG. 3: Pliase plane for the linearised evolution of d 
when A = — |. The line with positive gradient corre- 
sponds the growing mode for d shown in equation MUL 
and the line with negative gradient corresponds to the 
decaying mode (equation I4H '1 



To demonstrate more clearly the nature of the so- 
lutions to the linearised evolution equation for d, we 



present a phase-plane analysis. To transform equa- 
tion (|38|l to an autonomous form, we make a change 



of variable to s 



(3A+8) 



Int. Using this new time 



coordinate, the evolution equation becomes 
9 



S'^ 



3S' 



-A(A + 2)(5d = 0, 



(42) 



where a prime denotes the derivative with respect to 
s. If we put X — 5d and Y ^ S'^, then the system 
can be written 

X' = y , r' = ^A(A + 2)x-3r. (43) 

The system has only one equilibrium point, at X = 
Y = 0. The eigenvalues of the system around this 
point are ^ and "^('^+2) ^ -^jth the first negative 
and the second positive. Thus X = Y = is a. sad- 
dle point. The equilibrium lines of this saddle point 



are found to be 1" = ^X and Y 



-3(A+2) 



X. 



2 '""^ 2 

The equations (|43|l show that the phase trajectories 
go from left to right, i.e. towards the equilibrium 
line with positive gradient. Converting the equilib- 
rium lines back to our original coordinates confirms 
that they correspond to the solutions H40|l and (|41|l 
respectively, and that the growing mode l|in|l is an 
attractor line. The phase diagram for A = — | is 
shown in Figure |21 

Earlier we made the assumption that a and (j> dif- 
fered only by a small amount from the background 
solution. We now quantify this as 



2 h-i{Hat) 
(3A + 8) 



+ (5„(t), 4> = 



HHpt) 
V2b 



S4,it) (44) 



where (5^ and S,p are small relative to the back- 
ground. Now the linear reduction of the remaining 
evolution equations from |l7j is obtained as: 



3 So, (6A-f22)j^ 



12 



V2b t (3A 

-18\/2&(A + 4) 
(3A-K8)2 



Sa 



Sa 



8 Sr, 



ad.) 



1 Srf. 



(3A-t-: 
3(A- 



2) Sd, 



(3A- 



(45) 



6AV26 5^ 



(3A- 
0, 



t '3V2bt t (3A + 8)2 i2 



(46) 



2V2b Srj, 



t (3A 



(3A-k; 



t2 



Sri , 3(A + 2)<5d 



where 



ad*) 



t (3A + 8) ^2 
(1 - 



(47) 



{1_ (l_d^)(3A+12)/(3A+8} 
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The system (|2H|, 63), (EUl has a general 

solution of the form 

= rf^i-3(A+2)/(3A+8) ^ ^^^3A/(3A+8) ^ (48) 
5^ - 0lt-l + </.2t-6/(3A+8) 



0.49 



(3A + 8) 



3(A + 1) 
V26(3A + 4) 



(49) 



^3A/(3A+8) 



(3A- 



-6/(3A+8) 



6(A + 1) 



(3A + 8)(3A + 4) 



(50) 

^3A/(3A+8) ' 



where di, d2, 4>i, 4>2 and Cq are all arbitrary con- 
stants of integration. However, we can neglect as 
it is a gauge mode corresponding to a redefinition of 
the ordinary spatial coordinates. We can then see 
from the form of this solution that any perturbation 
to the background metric functions d, (j) will decay 
(in the linear regime at least). 
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FIG. 4: Numerical simulations of the non-linear radion 
evolution (dotted lines) are compared to the solution for 
the radion in the linear regime (solid line) 



B. Full non-linear low energy evolution 



Having completed the linear analysis, we solve the 
full non-linear system of low energy equations nu- 
merically, to see if the results above are substantially 
changed in the non-linear regime. The parameters 
chosen for this numerical analysis are 

(the value of 5 chosen for the previous numerical sim- 
ulation was too large to be valid in the low energy ap- 
proximations, so more appropriate values have been 
chosen). The background value of d is taken to be 
rf* = 0.5 and the corresponding tuned value of the 
potential on the second brane is a' ~ —0.023. Figure 
21 shows the numerical analysis (dotted curves) com- 
pared to the solution for the radion evolution from 
the linear regime, shown as the solid curve. We can 
see that at early time our linear approximation is 
reasonably accurate, but that as the collision is ap- 
proached the two curves are widely divergent. 

The dotted curves show the results for the nu- 
merical simulation with two different sets of initial 
conditions - one set of initial conditions correspond 
to the growing mode identified in the linear analy- 
sis, and hence to the brane trajectories identified in 
the exact solution, whereas the other set correspond 
to a mixture of the growing and decaying modes. 
The numerical solution for a and 4> under these ini- 
tial conditions shows that while they have identically 
the form of the background solution for the growing- 
mode conditions, the mixed-mode initial conditions 
set an initial perturbation away from the background 



solution which quickly decays to (or to the constant 
gauge mode in the case of a). This confirms the sta- 
bility of the background solution to perturbations 
of the brane position that was demonstrated in the 
linear regime. 



Accuracy of the low energy approximation 
and simulations 



Figures |S1 and show the proper distance between 
the two branes (given in equation (|33|l ) as obtained 
from the numerical integration of the low energy ef- 
fective equations and from the numerical solution 
for the brane trajectory in the exact bulk metric 
solution, plotted for several different values of the 
parameter 5. There is some difference between the 
exact and low energy curves, but we can see this 
difference decreasing as 5 decreases, indicating that 
the integration of the low energy equations is accu- 
rate in the limit 6^0. Figure El demonstrates the 
sensitivity of the system to even small changes in 5, 
whereas Figure El demonstrates the wide deviation 
of the approximate solution from the exact solution 
when 5 becomes larger. 

Figure [71 demonstrates the sensitivity of both so- 
lutions to a change in initial conditions. An increase 
in the initial perturbation in the distance between 
the branes dramatically decreases the time to the 
brane collision. 
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FIG. 5: The proper distance between the two branes 
as given by the exact solution (sohd hnes) and low- 
energy non- linear equations (dashed lines), plotted for 
S = 1/200 (lowest), S = 1/100 (middle) and 5 = 1/50 
(highest). 
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FIG. 6: The proper distance between the two branes as 
given by the exact solution (solid lines) and low-energy 
non- linear equations (dashed lines), plotted for 5=1/5 
(highest), 5 = 1/10 (middle) and S = 1/100 (lowest) 



IV. CONCLUSIONS 

In this paper we have demonstrated that for the 
bulk metric with scalar field given in 5] there is an 
instability in the position of a second brane. Solv- 
ing for the trajectory of this second brane shows that 
the instability can lead to a collision between the two 
branes. The scale factor and scalar field are shown 
to retain the same dependence on brane proper time 
regardless of the motion of the brane, and hence re- 
main regular at the collision. It was found in |^ 
that a general class of potentials for 2-brane systems 
with bulk scalar field will result in a singularity at 
any collision; however it was found in ^20], as here, 
that the exponential potential is exceptional and can 
tame the singularities at the collision. 

It is also worth noting that (for a countably in- 
finite subset of the possible values of the parame- 
ters b and A) the scalar field potential used here 
can be reproduced from a compactification of an 
empty higher dimensional space (18i il9]. with the di- 
mensionality determining the values of A and b. In 
this picture the collision is singular in that the bulk 
spacetime loses a dimension when the space between 
the brane disappears, however the induced metric of 
the brane and compactified space remains regular, 
and hence the effective scalar field derived from the 
compactified volume will remain regular also. 

In a complementary approach we considered the 
system using a low energy expansion developed in 
[l7|. This expansion reproduces the solution found 



in the previous section to good accuracy in the cor- 
rect limit (see Figure The low energy equations 
also identify a solution excluded by the methods 
used previously, in which the initial perturbation of 
the second brane necessitates a perturbation of the 
bulk metric away from its background. This solution 
quickly decays toward the original solution with the 
unperturbed bulk, confirming the conclusion of 0] 
that there can be no unstable homogeneous pertur- 
bations of the bulk metric. 

It would be interesting to see if the capability 
of solving the cosmological perturbations demon- 
strated in Koyama and Takahashi's one-brane solu- 
tion 6] is reproduced in this two-brane model. The 
perturbations would diverge as the branes approach 
collision, but it is possible to match divergent per- 
turbations across a collision in some scenarios (see 
e.g. HJl). If this is possible it would give an inter- 
esting parallel to the "Born-again braneworld" ^2Tj, 
the cyclic model (22| and other pre-big-bang scenar- 
ios. A microscopic particle-theoretic understanding 
of interactions at the collision would be helpful in 
determining the dynamics at the collision (whether 
the branes are likely pass through each other as in 
ll^, 112 , coalesce as in |2p| or disintegrate as in 
[23), as well as the behaviour of perturbations near 
the collision. However there are techniques available 
to evolve perturbations through brane collisions at 
a purely classical level |3 • 
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